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LONG-TERM RADIATION EFFECTS ON GaAs SOLAR CELL CHARACTERISTICS 

By ' . ■ ■ , ' 

J.H. Heinbockel‘ and M.J. Dbviak^ 

SUMMARY 

This report investigates preliminary design considerations 
which should be considered for a space experiment involving 
Gallium Arsenide (GaAs) solar cells. The electron radiation 
effects cn GaAs solar cells were conducted in a laboratory 
environment, and a statistical analysis of the data is presented. 

In order to augment the limited laboratory data, a theoretical 
investigation of the effect of radiation on GaiAs solar cells is 
also developed. The results of this study are empirical prediction 
equations which can be used to estimate the actual damage of 
electrical characteristics in a space environment. The e.xperimental 
and theoretical studies also indicate how GaAs solar cell parameters 
should be designed in order to withstand the effects of electron 
radiation damage. 


‘ Professor of Mathematics and Computing Sciences, School of Sciences 
and Health Professions, Old Dominion University, Norfolk, Va. 23503 

Assistant Professor of Mathematics and Computing Sciences, School 
of Sciences and Health Professions, Old Dominion University, Norfolk 
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INTRODUCTION AND STATEMENT OF THE PROBLEM 


This study was initiated to investigate preliminary design 
considerations for a space experiment on GaAs solar cells. The 
first steps in the design of a space experiment on GaAs solar 
cells are to identify the experiment goals and determine the need 
for a space experiment. The actual feasibility of a space experi- 
ment and the conceptual design of a flight package are not considered 
in this study. 

• 

In developing the goals of a, space experiment consideration 
should be given to the advantages which GaAs solar cells offer 
over other types of solar cells. These advantages are (1) higher 
AM-0 efficiency, (2) greater stability in a space radiation environ- 
ment and (3) operating capability up to 300® C. 

The need for a space experiment arose from two needs. First, 
laboratory measurements can be made which accurately compare 
different solar cell types. However, reliable data is needed to 
show how GaAs solar cells perform in a real space orbit rather than 
in an environment imperfectly simulated in a laboratory. Secondly, 
a functioning demonstration is needed to convince spacecraft designers 
that GaAs solar cells are ready for use in space power systems. 

The first experimental goal became the measurement of the power- 
generating capability of GaAs solar cells when exposed to a space 
radiation environment. Data on the effects of radiation on GaAs 
solar cells was just becoming available when this study was initiated. 
The data was generated through NASA/Langley Research Center and 

contractual efforts. Some ■ limited data was also available through 

independent programs . 

The second experimental goal was the measurement of power 
generation in GaAs solar cells at elevated temperature and various 
sun intensities in space. A review of the available data on elevated 
temperature operations of GaAs solar cells indicated that insufficient 
data existed to plan a 300®C space experiment. The need for such an 
experiment was demonstrated through intensive laboratory research. 
Because data is just now becoming available, the high temperature 
study is not considered in this report. 
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LIST OF SYMBOLS 

electron charge = 1.602192(10“^’) (Coulombl 

density of incident photons per unit time per unit 
bandwidth 

fraction of photons reflected from surface 
speed of light = 2.997 (10’) [ms“M 
fluence [electrons/cm^ ] 
number of acceptors on p-side 

number of donors on n-side 

intrinsic carrier density 

Boltzman's constant = 1.38(10’^’) [JK~*] 

Plane)? 's constant = 6.6262 (10“”') [Js] 

effective masses of holes and electrons [)cg] 

I 

electron rest mass = 9.109558(10“’*) [3?g] 

band gap energy [eV) 

temperature [*KJ 
permitivity of GaAs (F/cm] 

window thic)cness [cmj 

mobilities of minority carriers [cm^/volt-sec] 

“junction depth (ym] 

width of depletion region [cm] 
total cell thicl?ness [cm] 
recombination velocities [cm/sec] 

lifetimes of minority carriers [sec~*] 

diffusion coefficients [cm^/sec] 

,«.ceomG PASE «0T P.UMED ^ 


L , L , L diffusion lengths [cm) 
a g p 

a, 0 absorption coefficients for GaAs and GaAlAs [cm M 

J photocurrent density [mAmps/cm^] 

I short-circuit current density (mAmps/cm^ 1 

sc , 

SR spectral response 

P(X) average spectral irradiance at l.A.U. [W cm~^ ym“V] 

I 

E photon energy [eV] 

\ v/avelength [ym] 

open circuit voltage [V] 

theoretical maximvim open circuit voltage [V] 
y^ response of dependent variable 

0^ unknown parameters 

normally distributed random variable 


y, £ column vectors 

a standard deviation 


STATISTICAL ANAI.YSIS OF GaAs SOLAR CELL DATA 
Introduction 

The next several sections will discuss the effects of 1 MeV 

electron irradiation on the performing capabilities of GaAs solar 

cells. In the following analysis the cell junction depth Xi 

(2 levels) and fluence (4 levels) will be the independent 

variables. The ratios (I /I , ) of short-circuit 

sc SCO oc oco 

current and open-circuit voltage after irradiation to that before 
irradiation will be the measures of the . solar cell degradation 
and will be the dependent or response variables. The analysis 
consists of performing a multiple regression analysis to relate 




I /I and V /V to both Xi and j . The resulting 

sc SCO oc oco 

prediction equations then can be invoked together with estimates 
of electron flux in appropriate orbits to predict decreases in 
the short circuit-current and open-circuit voltages for various 
orbital periods. Residual analyses and tests for lack of fit 
are performed to check the adequacy of the models as well as 
underlying assumptions. 

Due to the complex functional relationships existing between 
the aforementioned variables, the models developed -are empirical 
in nature. Many modeJs were entertained before formulating the 
final models presented in this report. Many models were rejected 
due to a poor fit or intractable analysis. When used properly, 
the final models presented in this report will provide good pre- 
dictions for the response variables. These models fit the 
experimental data well and are also relatively concise. 

Description of Experiment 

The experimental units in the experiment consisted of 45 
heterofaced p-GaAlAs/p-GaAs/n-GaAs solar cells fabricated using 
the etch back epitaxy process. A detailed account of the growing 
process may be found in reference 1. Large area vacuum evaporated 
Sn-Aq contacts were used for electrical connection to the N-GaAs. 
The front finger contacts were spurted Pd-Ag. 

The solar cells were mounted on aluminum backing plates and 
then irradiated with 1 MeV electrons. Groups of cells were 
irradiated at fluences of 1C*’, lO*"*, 10*’ and 10*’ electrons/cm'. 
The short-circuit currents as well as the open-circuit voltages 
were measured before and after irradiation. 

Due to lose connections during the measuring process, some 
faulty data was obtained. Thirty-four measurements on short- 
circuit current and 32 measurements on open-circuit voltage 
remained after the faulty data was discarded. Each cell had a 
junction depth of either 1.5 or 4.0 microns and was subjected to 
one of the 4 fluences of 1 MeV electrons as mentioned above. 
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Experimental error explains different values of the dependent 
variable occurring at the same setting of the independent variables 
and occurs because of factors beyond the control of the experimenter. 
Some of these factors are the technique of growing the cells, 
heterogeneity of the experimental units, differences in the am.ount 
of irradiation administered, mechanical failures and measurement 
errors,. 

Results of Regression Analysis 


Let Y^, i=l,2,..., n denote a response or dependent variable 
that is measurable. We assume that a relationship of the following 
form exists between the independent variables and dependent variables 


= {Jo + 


+■ , i — 1,2,..., n 


where, 3o» ? i » . • • / are unknown parameters, 

known values of the k independent variables and s . is a normally 
distributed random variable with mean 0 and unknown variance 
[written ~ N(0,c^)]. 

The above model is a linear model, since it is linear in the 
parameters, and it can be represented in the matrix notation 


Y = X 3 + e 


( 1 ) 




where Y* = [Yi, Y:,..., Y^J , 3' = 13 


G t 


^ ~ [tl, £2, 


. . , £^ ] and 
n 


X is an n X (k + 1) matrix with typical element 


(x^^ = 1 for all i) 


^i] '710 " ' ■■■■ 1 

distributed normal variables with mean 


The £,*3 are independent and identically 


£ given by a‘I where I 


■ . The 
is the 


variance-covariance matrix of 
n X n identity matrix. 

. Assuming that X'X is nonsingular, the least squares estimator 
for 3 is given by (refs. 2,3): 


b = (X’X)"' X' Y 


( 2 ) 


= [Y'Y - b'X'Y] T [n - (k + 1]- ^ 

Although they will not be presented here, the esriniators b and 
possess many desirable properties from a statistical stand- 
point. The equation 



is called a prediction equation and may be used to predict future 
responses of the dependent variable at the desired setting 

of the independent variables. A (1 - a) 100 percent prediction 
interval for a future response Y at is given by 


where t. v is obtained from a table containing probabilities 
for students' -t distribution. 

In some instances the variance-covariance matrix of e is not 
o-I but instead o"V where V may have nonzero off diagonal 
elements (implying correlations between responses) or unequal^ 
diagonal elements (implying some responses are more variable than 
others). In these instances a weighted least square analysis is 
needed to determine equations (2), .(3) and (5), and these formulas 
must be accordingly modified to read 

b = (X’V-'X)-'X'V-‘Y 


2 = Y'V^ Y - b'X'V“‘ Y 


and 


1 , 

• /V til— I V \ ~ 1 V 1 


Y i V o[V. + X^MX'V-‘ X) X^J ' 


respectively, where 


V is Che variance of 

i 


(7) 


( 3 ) 


Y at 


For the short-circuit current data the following model is 
assumed: 

= 00 + + ^2^12 ^i ^ = 1,.../ 34 (9a) 

where = logio ^ = j' 


junction depth of cell 


and = logxb [<{)/10^®] with <j> the fluence. The matrix 

is assigned to be a diagonal matrix with diagonal elements 




and 


dll 

s= 

d2 2 

= 

0.0876 



ds 3 

= 

d4 4 

= 

d 1 3 / 1 3 

= 

1.413 

dl 4 r 1 ■» 

= 

• • • 

= 

1 9 / 1 9 

= 

0.3532 

d2 0 f 2 0 

= 

• • • 

= 

d2 3 / 2 3 

= 

0.1554 

d2 4 / 2 4 

= 

'* « « 

= 

d 3 4 / 3 4 

= 

1.413 


The weighted least squares estimate of 6 is computed to be 


and 0 = 0 . 138 


0.8S49 
b = I 0,1940 
1.075 

which gives a prediction equation of 


Y = 0.8649 + 0.194X1 + 1.075x2 


(9b) 


or 


sc 


SCO 


1 + /^lO®’®®*** 

\ Alo'v 


(9c) 


A plot of equations 9a to 9c is given in figure 1, 






Prediction Interval 

Let xi = 0.8 ym and = -2 (^ = 10‘‘»). Then a point predictor 
for a future response Y is given by -1.1299 and thus a point 
predxctor for is calculated to be 0.931. A 95 percent 

prediction interval for Y is given by 

-1.1299 ± 2.041 (0.138) [1.413 + X^(X’V“^X)-^X^]^/^ 

where 2.041 is the tabulated t-value, 1.413 is an estimate of V. 
and 0.138 is tf.e computed estimate of a. The above reduces to ^ 

-1.1299 ± 2.041(0.136) [1.413 + 0.112]*/^ 
or 

(-1.478, -0.782) 

which is the prediction interval for Y. This gives a 95 percent 
prediction interval for of (0.858, 0.967). 

The quantity Y^ - y , i = 1,2,. ..,34 is called the ith residual 
and gives the difference between the ith value of the response 
variable and the value of the prediction equation. Plotting the 
residuals against various quantities provides a graphical method 
of checking whether the prediction equation fits the data well and 
also whether the underlying assumptions for the analysis are violated 
Figure 2 gives the weighted residuals for the short-circuit current 
data. These residuals are plotted against Y, and, with the 
exception of one point, ^the residuals tend to fall in a horizontal 
bend around the Y-axis. ' If an unweighted analysis were performed, 
the variance of, the residuals would differ significantly and they 
would not have the horizontal band appearance which is desirable. 
Also, because there. is no curvature present in the residuals, the 
prediction equation fits the data well. 

The analysis of variance for the regression analysis is given 
in table 1. 
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Table 1. Analysis of variance for regression analysis of ^gc/j^sco' 
DF ■- degrees of freedom, SS = sura of squares of error. 


MS = 

mean square 

error, and F = 

F“statistic . 


Source 

DF 

SS 

MS 

F 

Regression 

2 

28.255 

14.128 

761.217 

Residuals 

31 

0.575 

0.019 


Total 

33 

28.830 

i 



In table 1 the F-statistic of 761.317 is highly significant, 
meaning that the junction depth and f luence terms should be included 
in the model. A more extensive analysis of variance, which decomposes 
the residual sum of squares into sim of squares pure error and sum 
of squares lack of fit, is given in table 2. 

Table 2. More extensive analysis of variance; DF = degrees of 

freedom, SS = sum of squares of error, MS = mean square 
error, and F = F-statistic. 


Source 

DF 

SS 

MS 

F 

Regression 

2 

28.255 

14.128 

761.317 

Residual 

31 

0.575 

0.019 


Lack of fit 

5 

0.071 

0.014 

0.74 

Pure error 

26 

0.504 

0.019 



The additional F-value in table 2 is a test statistic for 
testing for a lack of fit in the model. Had this value exceeded 
the tabulated value of Fo. os, 5 , 2 s = 2.59 then lack of fit could be 
concluded. However, since 0.74 < 2.59 there is no evidence of 
lack of fit. 


Open-Circuit Voltage Model 

A similar regression analysis was performed on the open-circuit 
voltage data. Several models were considered before a satisfactory 












r.y^.9 1 :■ I j 


fi.t to the data was found. The model of equation (9) with 

Y. = log 1 ,( 1 ^ 

\ oc 

residuals versus fluence suggested that a quadratic term in the 
fluence be added to the model. The resulting model is 


- 


was fitted to the data; however, a plot of 


« 3o + ^ “ 1 ,..., 32 ( 10 a) 

where = log i a ~ 1 ^* = junction depth of cell i 

and Xj^j = logi 9 ('j>/10' * 1 . Again it was necessary to use least 
squares with weights given by 


dll — ... = ds 6 “ 0.062, d77 « ... = dii,ii ~ 0.25 

dvi,xi = ... d\7,i7 = 1, di e , \ j = ... = d2x,2j * 25 

d2 2 , 2 2 “ ••• d 2 5 » 2 5 = 0.64, d 2 s , 2 s = ... ~ dx 9 , 2 9 = 1 

and daofsq =» ♦.. = d 3 2 »s 2 = 0.64. The least squares estimate 

of 3 was computed to be 


b * 


-0.6037 

0.098 

0.609 

0.060 


which gives a prediction equation of 


Y * -0.6037 + 0.098X1 + 0.609Xj + O.OGO-'ti 


( 10 b) 


The predicted values of j'^oc'^oco different junction 

depth-f luence combinations is illustrated in figure 3. An 
examination of the residuals in figure 4 shows no glaring 
abnormalities. 

The analysis of variance for the regression analysis is pre- 
sented in table 3. 
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irradiation at various junction depths. 
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Table 3. Analysis of variance for regression analysis of 


DF = 
MS = 

degrees of 
mean square 

freedom, SS = 
error, and F 

sum of squares of 
= F-statistic. 

error. 

Source 

DF 

SS 

MS 

F 

Regression 

3 

2.988 

0.996 

571.219 

Residuals 

28 

0.049 

0.002 


Lack of fit 

3 

0.007 

0.0023 

1.43 

Pure error 

25 

0.042 

0.0016 


Total 

31 

3.037 




In table 3 the F-statistic of 571.219 is good and the F-statistic 
testing for lavck of fit is not significant at a •= 0.05, which 
indicates there is no lack of fit. When the model was fitted 
without the quadratic term there was, at a = 0.05, evidence of lack 
of fit. 

SPACE EXPERIMENT CONSIDERATIONS 


The following is an outline of statistical considerations 
that roust be incorporated into an in- space experiment for comparing 
GaAs solar cells with other varieties of solar cells. In particular, 
we will discuss in a rather general way such things as: the design 

model and assumptions associated with an experiment; a discussion of 
the data analysis; and sample size considerations. Finally we will 
summarize the proposed techniques. 

Let ^sc^^sco' denotes the ratio of the short-circuit current 

after irradiation to that before irradiation, be a measure of the 
solar cell degradation. (We could also use maximum power ratio or 
open-circuit voltage ratio as a measure of solar cell degradation.) 

Let I denote the short-circuit current measured at some point 

S w 

during or at the conclusion of a space flight. We let p denote 
the number of types of cells included in the experiment. The cells 
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will be either GaAs or and some may be chosen to be annealed. 

In any case there will be p different types of solar cells under 
consideration. 

Notationally we let i = 1,2,..., p be the average value 

of ^sc^^sco cell type i at the designated orbit where the 

average is based on a hypothetical population of cells. Actual 
responses of '^c'^^sco differ considerably from the corresponding 

due to experimental error which is composed of: measurement 

error, cell heterogeneity, faulty connections, and different amounts 
or types of irradiation. Even though experimental errors exist we 
can still meike inferences about the v^'s and more specifically 
about - Pj, i j. 

Previous sections are of a rather general nature, partially 
because little is known concerning the proposed flight at this 
time. More specific recommendations could be made as the decisions 
on the proposed flight become finalized. The following decisions 
need to be made concerning an experimental design: 

I 

(a) What area is available for solar cells onboard the space 
vehicle and the number of cells that can be mounted on this area, and 

(b) whether the experiment is to be a multi-factor experiment 

(with three factors: (1) junction depth (2) type of cell, or 

GaAs, and (3) whether cells are annealed or not annealed) . 

If a multi-factor experiment is desired the considerations 
discussed earlier should be changed accordingly. However, the 
previous discussion on sample size will still be applicable. In 
any case, the goals of a proposed experiment must be specifically 
stated prior to the experimental design stage. 

Experimental Design 

The only apparent source of variation , other than those listed 
above, contributing to experimental error is due to solar cell 
(trea'jnent) differences. Therefore a completely randomized design 
is a logical choice. This means that if n cells (assuming that 
the same number of cells of each type will be included) should 


be randomly chosen from all available cells of each type, there 
will result np experimental units in the experiment. A table 
of random numbers (or similar device) should be consulted to 
guarantee random samples. 

The probabilistic model associated with a completely randomized 
design is given by 

y^j ~ + ^ij ^ ~ 1/2,. ../P j — 1,2,..., n 

where is the jth response (value of ^gc^^sco some other 

measure) at treatment i, the y^'.s are the unknown treatment 
means, and is the random error associated with the measuremen 

Y^j. It is assumed that for all i and j, is a normal 

random variable with mean and variance o and ct^, respectively. 
It is also assvimed that any pair of e's are uncorrelated 
random variables. 

If it is suspected that the variance of the response changes 
with treatments, then an appropriate transformation of the data 
should be made in order to stabilize the variance. This technique 
was used in the regression analysis of the terrestial data. If 
the data is transformed it is noted that the resulting confidence 
intervals apply to the transformed and not to the original 
parameters. 

ANALYSIS OF DATA 

/ 

The first step in the statistical analysis is to perform a 
one-way analysis of variance on the data. This entails decomposing 
the total sum of squares of error | total SS = (Y. . - Y) ^ j into 

7 ' - - ^ 

treatment sums of squares I SST = 2n(Y^ - Y) 1 and error sum of 

squares (SSE = total SS minus treatment sum of squares) , where 

Y = 22 Y. .in p and Y. is the mean of the observations for 
ij 

cell-type i. The results are summarized in the following analysis 
of variance (ANOVA) table; 



ANOVA 


Source DF 

Treatments p-1 

Error np-p 

Total np-1 


SS MS 

SST MST 

SSE I*!SE 

Total SS 


F 


MST/MSE 


The abbreviation DF appearing in the table denotes "degrees of 
freedom," the parameter of a random variable. The MS column 

is obtained by dividing SS by tlie appropriate DF. If F = MST/MSE 
is "small" (smaller than a te±iulated F value, based on p-1 and 
V - np-p degrees of freedom) then it is concluded that no difference 
exists in the treatment means and the analysis is complete. If 
the F-statistic is larger than the tabulated value it is concluded 
that not all the y^'s are equal. 

When the F-statistic is significant (large) it is usually 
desirable to determine just where the differences lie in the Uj^’s. 
The statistical techniques appropriate for this type of analysis 
are collectively )cnown as multiple comparison procedures (see ref. 

4) . Some of the multiple comparison methods in common use include: 
the Least Significant Difference, Scheffe's Method, Tukey’s Method, 
the Student-Newman-Keuls Method, Duncan's Method, and Dunnett's 
Method. The appropriate method in any situation depends upon the 
overall goals of the experiment and should be chosen accordingly. 
Although it is not obvious (see ref. 5) which procedure is appro- 
priate in some experiments, the TuJcey procedure is a li)cely choice 
for the experiment under consideration. It. ranics high in controlling 
the experimentwise error fate (ref. 5) and is particularly good 
if we wish find confidence intervals for the differences in 
treatment means. Briefly, the procedure may be outlined as follows: 
set an appropriate confidence coefficient 1 - a. (This is the 
probability that all the confidence intervals constructed will 
contain the respective parameters.) Obtain the quantity Q(1 - a, 
p, v) from special studentized range tables (ref. 5) . Compute 
Q(1 - a, p, v) (MSE/n) = D. Construct a confidence interval of 
any difference by calculating 
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LJ 


Y. - Y. 

> 3 


D. 


Do this for any dif_arences of interest. If a particular interval 
does not contain the point "zero," the means are concluded to be 
different. 

Sample Size Considerations 

Choosing the optimum number of cells of each type to include 
in the experiment is an important consideration in -the design of 
the experiment. If cost, weight, and available space restrictions 
did not exist, the more cells included in the experiment the 
better the inferences that, could be made. Since this is not the 
case, the sample sizes should be chosen to satisfy some criterion 
(or criteria) specified by the experimenter. One such criterion 
is to choose the n^'s so that the confidence interval widths 
(or half -widths) for the differences in the means is less than 
some bound (specified by the experimenter) with probability 1 - a. 
The choice of the appropriate bound is important. If chosen to be 
too large, the confidence intervals will be too wide and imprecise 

I* 

inferences will result. Too small a bound will yield inordinately 
large sample sizes. Some compromise is needed. For example, it 
may be desirable that the half-width of any confidence interval 
for ^ Wj not exceed 0.03 and that this will occur with 


probability 1 - a = 0.90. The optimum value for n 
obtained by solving 


can be 


0.03 = Q(0.90, p, V) (MSE/n) 


V. 


for n. For illustrative purposes assume that p = 8 and 

= 10^" (or some equivalent value when considering other types 
of irradiation). Also assume a junction depth of 1.5. Then 


approximating 

gives 


MSE 


V; 


bv 0.022 (based on the terrestial data) 


0.03 = Q(0.90, 8, 3 (n - 1) ) 


0.022 


/ 

/■ 


\ 




Consulting a studentized range table and solving for n gives 
the optiir.al n as approximately 8. Therefore if 8 cells of each 
type are included in the experiment the probability is 0.90 that 
all half-widths of the confidence intervals will not exceed 0.03. 

Actually, during the space flight a cell may incur any of a 
number (say k) of types of malfuncticns which may render the 
cell inoperable or make it impossible to obtain a measurement on 
the cell. In this case it would be reasonable to include n' cells 
of each type in the experiment where n* is some n^imber larger than 
the optimal number derived in the previous paragraphs. A method 
for approximating n' is now given. 

Notationally, let p^ i = 1, 2 ,..., k denote the probability 
that any cell incurs malfunction i during flight. Assuming 
that the types of malfunctions are independent events, the proba- 
bility that a cell incurs no malfunctions during flight is given 
by 

q = -£-(1 - Pi > . 

i=l 

Let 

I I if cell i has a defect free flight 
0 otherwise . 

Assuming the 's are mutually independent random variables, 
n' 

T = ^ = number of cells operating satisfactorily at the end 

of the flight,, and is a binomial random variable with 

parameters n' and q. A conservative rule for calculating n’ 
would be to find the binomial distribution whose probability of 
being at least n (optimal value determined previously) is 0.90. 
This would provide assurance that a sufficient number of cells 
will remain operable. The Pj^'s can be estimated if results 
from similar experiments are available. If not, the experimenter 
may have to subjectively estimate the p^'s. At least conservative 
upper bounds for the p's can certainly be given. 





Electron Irradiation In Space 


Although the data from radiation effects on GaAs solar cells 
is limited, we will show how to predict degradation in a space 
environment due to electron irradiation. It must be emphasized 
that the cells tested in the labcvatory were not radiation hardened 
and so the predictions made here will be underestimates for. life 
expectancy in an actual space flight. The prediction equations 
are based upon deep junction cells which are expected to degrade 
rapidly. 

From the laboratory data of solar cell degradation versus 
fluence analyzed earlier, we have obtained the prediction equation 


'te - ‘) ■ 


■J.8649 + 0.194xj + 1.075 log ((J>/10^ ® ) 


where <(i is the fluence [1 MeV e/cm^] and x^ is the junction 
depth [ym] . For a 25 percent reduction in the short-circuit 


current 


^sc^^sco “ 0*75, equation (11) , can be expressed 


'<}) = " 0 « 180 SXj 


We assume that = Qt where Q = Q(E > 0.5) is the integrated 
flux of electrons with energies greater than 0.5 MeV, which has 
the dimensions [electrons/cm^-sec] . Values of Q are obtained 
from figure 5, (from ref. 7) for various orbits. For synchronous 
orbits the altitude is between L = 5 and L = 7 earth radii and the 
flux Q varies between 5(10®) and 2(10®) electrons/cm’-sec. 

Using the equation (12) with x^ = 0.5lym], table 4 gives estimates 
for the predicted time t (days) for a solar cell to degrade by 25 
percent. 

















Table 4. Time (days) for GaAs solar cell to degrade 
by 25 percent for a junction depth of 0.5. 


Earth Radii 


Electron Flux with Energies 
Above 0.5 MeV 


Times (Days ) to 
Degrade 25% 


L 

5 

6 
7 


Q 

8 ( 10 ®) 

5(1Q«) 

2 ( 10 ®) 


^ 

663 

1061 

2653 


The values of t in table 4 are lower bounds on the expected 
time to degrade 25 percent and are based upon the prediction 
equation obtained from laboratory data on deep junction cells. This 
prediction equation does not take into account the temperature 
effects on the operating characteristics of the GaAs solar cell. 

GaAs solar cells can operate up to 300®C, and laboratory data is 
just now being obtained (ref. 8) which demonstrates that radiation 
damage caused by electrons can be removed with annealing. Thus, 
the predicted times to degrade by 25 percent are probably way too 
low. It is expected that with proper annealirlg of radiation 
damage encountered in an orbit there can be imposed a significant 
reduction in the expected radiation damage. However, further 
research into this important area is needed to verify this 
conjecture. 

Summary 

The completely randomized design previously discussed (see 
"Experimental Design") was chosen since only one source of variability 
(difference in treatments) was obvious. A one-way analysis of 
variance is then the appropriate first step in the analysis. If 
the experimenter so desires, the experiment could be considered as 
a two-factor experiment v/ith factors (1) type of cell and (2) 
amount of annealing. A two-way analysis of variance would then be 
the first step in the analysis (ref. 6) . The sample size dis- 
cussion would be similar, however. 


THEORETICAL CONSIDERATIONS 

In order to augmenc the limited data on radiation damage to 
GaAs solar cells, a modeling approach is presented in order to com- 
pare theory and exper j-mental results. This modeling approach con- 
sists of developing a mathematical model based upon the diffusion 
of holes and electrons within the various regions of the GaAs solar 
cell. The resulting model contains various physical parameters such 
as diffusion lengths, mobilities, recombination velocities, etc. It 
is assumed that the degradation in the electrical properties of a 
solar cell that has been exposed to an electron irradiation can be 
explained by determining the effect of the radiation on these para- 
meters (i.e. degradation of device parameters) - 

Background Material 

The following is a summary of existing models of semiconductor 
behavior (refs. 9 to 17). Consider the separate regions of p- and 
n-type semiconductor materials which are combined to form a junction 
as in figure 6. The n-type material has a large density of electrons, 
and the p-type material has a large density of holes. In the n-type 
material the electrons are referred to as majority carriers, and any 
holes created in the n-type material are referred to as minority 
carriers, in the p-type material the holes art the majority carriers 
and any electrons that are freed from bonds are called minority 

When the n- and p-type materials are combined there is a 
diffusion of holes from the p-region to the n-recion, and simul- 
taneously there is a diffusion of electrons from the n- to the p- 
region. This diffusion process continues until an equilibrium 
state has been achieved. As the electrons diffuse from n to p 
they leave behind donor ions nJ in the n material, and holes 
diffusing from p to . n leave behind acceptor ions N~. Thus atoms 
on either side of the junction become ionized. The charge associ- 
ated with an ionized atom is not free to move about within the crystal 
structure, and thus there is produced an electric field E which is 
directed from the positive charge toward the negative charge. This 
electric field E appears in some region W about the junction which 
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p-type 



n-type 


®, 0= Bound charges of impurity atoms in n- 
and p-type material 

+ , - = Positive (holes) , negative (electrons) 




= Conduction energy 


E = Valence energy 


E^ = Fermi level 


Figure 6. (a) Depletion region, (b) electrostatic potential, 

and (c) energy band diagram of p- or n-type material, 



is called the depletion or space charge region, and there is 
an equilibrium potential difference “ '^n ~ '^p region 

W, where v is the electrostatic potential and “ 

At equilibrium the drift and diffusion components of the hole 
current become balanced. Thus, letting J denote the current 
density, we have at equilibrium 

■^drift "^diffusion ~ ® 

q Up p(x)E - q Dp if = 0 (13) 

where p = p(x) = hole density tcm“^) 

Dp = diffusion coefficient [cm^/sec] 

Up = hole mobility [cm^/ (volt-sec) ] n 

q = electron charge tcoulomb] 


Jp(X) = 
Jp(X) - 


The diffusion coefficient and mobility are related by the 
Einstein relation 

Dp = Up KT/q . (14) 


where k = Boltzman's constant. (J/®K] and T = temperature [®K] . 
From the equations (13) and (14) we obtain the integral 



dv 

dx 


dx 



1 ^ 
p dx 


dx 


which gives 


za 

KT 




( 15 ) 


where p , p 
(electrons) 


, (n n ) are the equilibrium densities of holes 

P n, p 

on either side of the junction. At equilibrium 
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= nt = p Nj 
an 1 a 


where is called the intrinsic carrier density, if there 

ara acceptors on the p-side and donors on the n-side 

we can write equation (15) as 


Vo - f ‘n 


nJ/N. 

^ 1 '■ 


= £2 « 
" V"L 


This is/ the theoretical maximum open-circuit voltage that can 
be achieved. From equations (15) and (16) we let 




and n 




where V^' is the theoretical maximum electrostatic potential 

difference across the depletion region and P denotes the value 

, ''o 

P„ at X = “W and h denotes the value of n at x = -W . 
" P Pq P n 

We th^n have in terms of the actual potential difference V - V 

r . C 


across the junction that 


-q(V^ - V) 


qV/KT 

P^ e ^ and 
o 


-q(V^ - V) 


n = n e 
P n 




qV/KT 


The carrier density varies with temperature and band gap energv 
according to the relation 


n. = ‘ (.M*M*) ' 

i V j ^2 y n p' 


A -Eg/2KT 
e ’ 


•where M*, M* are the effective masses of holes and electrons 

in the semiconductor material. For M* ~ M* M , M* = M* M , M* 

n e o p hoe 

0.068, = 0.53 and = the mass of an electron, equation 

(19) can be written as 
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n? = (M*M*) (23.3172) (10*® 

1 e n 


- 0.84(10*®)T*e“%^^'^ 


where the band gap energy E^, in electron volts, can be approxi- 
mated by 

E = 1 S22 - ^•8g0~-)'r^ 
g T + 300 

E = 1.43 ev @ T = 300»K (21) 

9 

where T = temperature [®K] , K = Boltzman's constant [eV/®K] . 

The width W of the depletion region is W = and 

can be determined by examining the charge density within the 
depletion region. t equilibrium we have 

qAWpNa = qAW^N^ (22) 

where A is the cross-sectional area of the solar cell. Using 
Gauss's law, the gradient of the electric field can be related 
to the space charge at any point x about the center of the 
space charge region and we can write (see fig. 7) : 


ii - f »d' 0 i i «„ 


II = ^ N , - W < X < 0 

ax e a p — — 

where e = ^r^o permitivity of the material and 

= 8.85(10"*“) F/cm is the permitivity of free space (for 
GaAs = 10.9). 

Integrating the equations (23) gives 


E = ^ N ,W = — N W 
o e d n G a p 





ELECTRONS 


(b) 

Pp 

n 

n 

Majority carrier concentration 
(cm”’] 


Pn 

n 

P 

Minority carrier concentration 
(cm"’ ] 


n 

Po 

n 

Po 

Minority carrier concentration 
at thermal equilibrium (cm"’] 

n- type 


p-type 

Impurity concentration (cm“’] 


acceptor 


Figure 1 . (a) Deplet’^cn region and (b) symbolism. 
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Now, -V^ =/ " E(x)dx = I E^W = 


-w. 


“if N ,W W 
2 e an 


where W = W + W and from equation (24) N,W 
“ P an 

equations can be combined to give 


= 


2eV 


\ ^a^d 7 


N W SO that these 
a p 


(25) 


In the case of nonequilibrium where there is a voltage V applied, 
the altered value of the electrostatic potential is - V and 
equation (25) can be written 


a _ 


= 


V / 


(26) 


Gaj^_^ Solar Cell 

With reference to figxire 8, a model of a Gaj^_^ A1 As solar 
cell, the various processes of interest are that: 

1. Some photons are reflected from the surface. 

2- Photons with short wavelengths Igenerate hole-electron 
pairs close to the surface of the cell. Some of these pairs 
quickly recombine, and others are injected into the p-GaAs region. 


3. Photons generate hole-electron pairs in the p-GaAs. 

Here electrons move toward the depletion region and contribute 
to the short-circuit current at the junction edge. 

4. Holes generated in the n-GaAs region diffuse toward 
junction and contribute to the short— circuit current. 

Proceeding. as in references 9 and 10, the following equations 
for the photocurrent can be developed, when incident photons 
have an energy greater than the band gap energy, absorption of 
the photons causes electrons to be raised in energy from the 
valence band to the conduction band E^ creating hole- 

electron pairs. We will be concerned with excess minority 
carriers (holes on the n-side and electrons of the p-side) which 
diffuse to the edges of the depletion region before thev recombine. 
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It is these minority carriers which produce a photocurrent, 
photovoltage and consequently solar cell power. Let p^ - p^ 

denote the excess hole density in the n-GaAs material and let 

n_ - n denote the excess electron density in the p-GaAs 
^ ^o 

material. At a surface, the minority carrier current density 
is given by Jg^^^face' '^surface = - p^^^ n-type 

material [amp/cm^l = qS^^n^ ~ n^ p-type material 


with q ~ electron (hole) charge [coulombs] 

S = surface recombination velocity of holes 
Icm/sec] 

S = surface recombination velocity of electrons 
” [cn\/sec] 

p^ - = hole density (cm“M 

n - n = electron density [cm“'*] 

P Pq 

Consider the hole flux per unit volume entering and leaving a 
volume element of the solar cell. We have 


3 / Urn rJ(x) - J(x -^'Ax) T 

atV^n “ ^nj~ -*■ oL qAx J 


g + G 
^r r 


(rate of 
hole 

build-up) 


(increase 
in hole 
density 
per unit 
time) 


(recombination -t- (generation rate) 
rate) 


Where J is the current density of the holes. 
Continuity equation for holes as 


This gives the 


) , dj 

t q dx r 


Similarly, the continuity equation for electrons is 


3n , dJ 

P _ rJ; n _ ^ ^ 

3t ~ q dx ^r ^r 
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Let 3, a denote the absorption coefficients for the GaAlAs 
and GaAs regions respectively, then the generation rate for hole- 
electron pairs can be expressed as 


- R) e 0 < X < D 

a O — — 


Gg = a^oU - R) e”®°e"^ , D < x < d + 


Gp = ciF^d - R) , d + x^‘ < x < h (29) 


where R = fraction of light reflected from the surface of the solar 
cell and = number of photons incident per square centimeter. 

The recombination rates in the various regions are expressed as 


5a = 




n n 
P - P, 


, T. = lifetime of an electron in GaAlAs 

a ' 


n \ . 

P “ Po l 

)• “ lifetime of an electron in p-GaAs 

9 ' ^ . . I.'. 


^P 


P P 
‘^n - *^n. 


1, T = lifetime of a hole in n-GaAs 

P / ^ 


For steady-state conditions, the continuity equations for the 
different regions become: 

For electrons in the p- type GaAlAs material 


q dx a V / 


Ola) 


For electrons in the p-type , GaAs. material 





For holes in the n-type GaAs material 
, dJ 


- Pri 


- -r^ + G - 
q dx p 


= 0 


(31c) 


’ ^ T 

where - — -r- = divergence of current = increase in electron 
q dx 

(ho*e) density per unit volume per unit time, and G = generation 
rate. 

* 

The hole and electron diffusion currents in these regions are 


Ja = - q D 


a dx 


J = - q 

g ^ 


/n ''i 

p - Poj 


_ d (n n \ 
g dx \ P Pc/ 


J = q D ~ 
p ^ p dx 


- \) 


(32a) 


(32b) 


(32c) 


»{ 

1'^ 




where D , D / D are the diffusion coefficients for the GaAlAs, 
ago 

P-GclAs and n-GaAs regions, denoted by che subscripts a, g and 

P- 

From equations (31a) and (32a) we obtain the diffusion equation 
for the GaAlAs region as 


(n n) (n nj 

P - Po/ \ P - Po/ . , 


+ »(1 - R)F e " = 0,0^x<_D 


which is subject CO the boundary conditions 

qs/n - n \ = qD - n \ at :■ 

P a d-x ^ p pj 


s ■ 0 


ana 


n -n=0 atx=D. 
P 


(33) 
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From the solution of equation (37) the photbcurrent at 
X = D + Xj is given by 


•^D+Xj = - qo ^ 


g Sf("p - "pjl* = D + 


X . 

J 


Solving equation (37) cind calculating equation (38) gives 

J, 


'D+x. 


-qF^d - R)aL e~^^ Tfz 

+ C.L e-“^j 

a^L‘ - , 1^3 g 




where 


f 2 =? e 


-ax. . /^] fq\ 


I - ^ 


fa ~ — ^ ^ sinh(r-^l + coshlr-^ 

g \ g/ \^g> 


and 


Lg = yOgTg = diffusion length in p-GaAs 

For the base region, equations (31c) and (32c) give the 
diffusion equation 

D + „F (1 . R)e-eV“’'je-“(>^ 

^ dx D ° 


= 0, D + Xj +W£X<H 


which is svbject to the boundary conditions 


- p^ =0 at X = D + X + w 

o ^ 


V(Pn - PnJ - -=p ^(^n-%) ^ 


= H. 


( 33 ) 

(39) 

(40) ' 

(41) 

(D+x^ )) 


(42) 
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Solution of equation (42) gives the photocurrent at x 
as 


Jo + X. + w . qDp 

We calculate this current to be 


+W 


D + Xj + W 


(43) 


D+x .+W 
3 


qF (1 - R)oiL 

9 E 

- 1 


[■ . - f -] 

t? 


(44) 


where 


= 


s inh^ 

Sj) . -v- 

T,* S L 

•aH* p p 

D 

P 

- --'I 





II 

, D T and 

1 P P 

H* = H - 

(D + Xj + W) . 


(45) 


For large S we have fi, = 1 and equation (45) reduce'' 

to P 


Jd + Xj + W = 


qF^d - R)aLp e~^^e~°^3 e~°^ 

cLp + 1 


(46) 


Equations (35) and (39) have removable singularities where 
= 1 and = 1 ■ respectively. Introducing the functions 


f(x,y) = 


-y/2 , X = 1 

^ e^sinh(u), x ^ I, u = ^ ^ 


(47) 


and 




-qF^d - R)X 


= Y s ' inh 2 V cosh 2 |<y [fU,z) + 


X + 1 L 


11 - Y) 


- ll) 


( 48 ) 


Equations (35) and (39) can be represented in the alternate forms 



Some photocurrent collection also occurs in the depletion 
region. If the electric field in this region is high enough, 
then tlie minority carriers will be accelerated out of the region 
before they can recombine. The photocufrent contribution from this 
type of excitation is simply the number of photons absorbed and is 
given by ’ 


= photon density entering - photon density leaving 
= qF^d - R)e"®°e'°‘^3 [1 - 3 "“^] (51) 


The total photocurrent is the sum of the photocurrents given by 
equations (46), (50), and (51) and is denoted by J^. We have 

=>D+X.<*> ^ Vx +W<^> <52) 

Spectral Response-Short-Circuit Current 


The spectral response of the solar cell is a function of wave- 
length and is given by 



where in equation (52) the absorption coefficients a and 8 
are fionctions of wavelength A [microns] and p the percentage 
of Al in Gaj^_pAlpAs. These absorption coefficients are 
approximated by the relations : 


5.5074 (10-) VT 
1.2824 (10**) (E 


, E £ 1.38 1 

1.3779/E, 1.28 < E < 1.50 
0.5) , 1.50 < E 


0 E < El = 1.38 + 0.9535P 
6 = f I (E) , El < E £ Ea = 1.5 + 1.4 53Sp 
fa (E) , Ea < E 


where E = 1.2402/A is the photon energy [eV] at wavelength 
A [microns] and 


/ (6 ) (p) (Ea)\/ \ 

fa(E) = ^1 j\^1.2824 (10'*y(E - 0.5)^** 


(56a) 


fi(E) = 10 exp{Ki (E - El)} + 523p sin 


[■(E - Ej)tt1 

L Ea - El J 


(56b) 


with 


1 , /f2(Ea)\ 

" lEa - El) 


The absorption coefficients a and 3 are illustrated in 

figure (9) for various values of the percentage of aluminum p. 

The parameter S, raises (S, . near 1) and lowers (3,. near 0.25) 
iim xim iim 

the GaAlAs absorption curve- For this study = 0.75. 









The photocurrent at each wavelength A fmicrons] is given by 


Jp(A) = qF^d - R)SR(A) 


( 57 ) 


and the total photocurrent contributed from all wavelengths 
greater than or equal to Xi is the continuous sxm given by 


sc 


= q y “ R)SR(A)dA 


(58) 


which may also be expressed in the form 


sc 


= - 3 _ 
he 


L 


(1 - R)p(A) ‘A* SR(A)dA 


(59) 


where h = Planck's constant, c = speed of light, and p(A) 
tWcm“^um“^] is the average solar spectral irradiance at l.A.U. 
which has a solar constant of 135. 3 [mwcm“^ J . Values of p(A) 
can be found in references 10 and 11. 

The magnitude of the short-circuit current density at wave- 
length A is given by equation (58) and is 


■^sc^^^ “ '^p ~ “^D+Xj+W 

If there were no losses by reflection or recombination the 
maximvim current density would be 


(60) 


J , = qF 

sc max ^ o 


(61) 


Hence, the collection efficiency or quantum efficiency Q(A) 
is given by 


Q(X) = 


sc 


sc max 


■^D+x . + '^D+x . + W + J„ 
3 3 W 


(62) 


The short-circuit current density can be expressed in terms 
of the quantum efficiency by the relation 


li 


■'so “ “J ? F^Q(X)dX 
‘^1 


(63) 


Injected Currents and Recombination Currents 
For the dark current associated with the solar cell we have: 
(i) Electrons from the n-side injected to the p-side 

1 ("P ’ '"Po) 


q dx 


= 0 


(ii) Holes from the p-side injected to. the n-side 
q dx T 


We also have the current equations 

“^n dx(”p " 

^^p dx ^Pn 

Equations (64) and (66) produce the di.r . asion equation 
d' , ("p - "Po) 

D /n - \ — = 0 , D < X < 

^ dx^( P Pq) ‘'^g 

which is subject to the boundary conditions 


D + X . 
3 


s3V 


"p = p® pp 


D + X . 
3 


and 


j_/n - n \ = D :^/n - n \ 

9\P Pq) g dx\ p 


at X = D 


Equations (65) and (67) produce the diffusion equation 


(64) 


(65) 


( 66 ) 

(67) 


( 68 ) 
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xi 

:.a 

^ •• 


-.I, 


• . 1 ' 


i / 










“P dx"(^" ' %) 


(Pn - ^Rq) 

i -2Z = 0 , D + x. + W£x<H 


which is subject to the boundary conditions 


p =s p e at X = D + X . + W 


®p(^n ■ PnJ “ "°p dx(Pn " PnJ 


at X = H 


From equations (68) and (69) the injection current is given 


■^inj “ dx (Pn " Pn^j U =D+5c ^+W ^ dx^^p ' lx=D-hc ^ 

Solving equations (68) and (69) produces the injection current 

j. . = j [e - 1 ] -. (■ 

m3 o 


where 


( qD n^ , S L \ /qD n? H 

^ = W' - 


(D + x^ 
L_ 


W) S L 


-^1 ,7 = , 


with G = G(5,n,T) given by 

rtc „ ^ N _ c rsinh(n) ;+ t cosh(n)~ ] 
oK»n» ) - sinh(n) + cosh(n)J 


The recombination current in the depletion region has the 
maximum value given by (ref. 10) : 

im^WKT sinh(^^) 




(V^ - V)_ 




( 74 ) 


where is given by equation (17) and W is obtained from 

equation (26) . 

The total current is then given by 


J = - J. . - j 

sc inj rec 


( 75 ) 


The open-circuit voltage is obtained from aquations (59), (71), 
(74) and (75) by setting J = 0 in equation (75) and solving 

for the voltage V = ThtJ power output of the solar cell is 

given by 


p = JV 


(76) 


where J is obtained from equation (75) and at maximum power 

^ = '^max' ^ = '^max' = Pmax* ' 

be obtained by solving the equation 


P “ Pmax* voltage at maximum power can 


If = ^ * V § = 0 


(77) 


for the voltage, V = v 


max' 


The efficiency of the solar cell is given by 


EFF = n = 

Pin^t 


P A 
^max a 

(135. 3)A, 


(78) 


where is the total area of the solar cell and A is the 
actual area which equals the total area minus the finger area. 
The fill factor is calculated from the relation 


FF = 


max max 

J V 
sc oc 


( 79 ) 



ANALYSIS OF DEVICE EQUATIONS AND COMPARISON WITH EXPERIMENTAL DATA 


A computer program was written which incorporates the equations 

describing the solar cell behavior. In this progrcUTi the spectral 

response as determined by equation (53) was treated as a function 

of variables L , S , L , S , y and L . In addition, the 
a a g g n p 

mobility y^, y^ of the minority carriers in the p-region and 
n-region are related to the diffusion constants of these regions by 
the Einstein relations 



I ^ and = u ^ 
n q P P q 


(SO) 


The diffusion lengths L. , L„, L , the diffusion constants D;., 

“ 9 P ® 

D„, D , and lifetimes t,, i in the respective regions 

g p agp 

GaAlAs, p-GaAs, and n-GaAs are related by 


L? = D^T^ (i = a,g or p) 


(81) 


The "other" parameters occurring in the equations which model 
the solar cell response were assigned the nominal values which 
are given in table 5. 


The measured values of the spectral responses for each cell 
were normalized and compared with normalized values of the predicted 
spectral responses. Measured values of the open-circuit voltage 
and short-circuit current were also compared with the predicted 
values from the model. In particular, a Marguardt optimization 
algorithm was employed and the computer was asked to choose 
values of the parameters L^, S^, Sg, y^, which would 

minimize the sums of squares errors as well as minimizing the 
absolute values of the following percent errors: 


fi(Xi) 


Error (i) = 
(100) , i = 1, 


SR(X^) measured - SR(X^) predicted 
SR(X^) predicted 


(82) 
(cont ’d) 
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Table 5 


Nominal, values assigned to GaAs solar cell 
parameters (before irradiation values) . 


L = 2.4(10 cm 

Si 




Lg = 3(10“'*) cm 




L = 3(10*") cm 
P 




N^ = ( 10 ^®) 

a 


N^ = 1.3(10^^) 


D = 0.4 (10"“) cm 


Xj = l.SdO’”*) cm and 4.0(10““) cm 


= 10® cm/sec 


S = i0“ cm/sec 

g 


p = 0.83% aluminum 


Up = 300 


= 3300 i 


T = 300 


T = 1.0(10"®) 








f 1 0 (x^) = Error (10) 
f 11 (x_j^) = Error (11) 



where A^(i = 1,...,9) were the wavelengths 0.45 to 0.85 in steps 

of 0.05 rr.icrons. In the computer program, the variables x^(i = 1, 

..., 6) = (L ,S ,S ,u ,L ) were restricted to lie within certain 
& ^ n p . 

regions ^imax* accomplished *by introducing 

the dummy variables u. (i = 1 ,..., 6) which were related to x. by 
the relations 


( X‘ X \ ■ / X X \ 

:max + imin ] . I imax - imin \ 

~“ T “ / \ ~~2 / 


tanh (u^ ) 


The variables u^^(i = 1,..., 6) were then allowed to vary over the 
range < u^ < «>. 

The Marguardt procedure utilized can be represented by the vec- 
tor equation 

= {x^} - |^j'^(x^)J(x^) + j'^(x^)F(x^) 

where 


F(x*^) = {(Col(f 1 (x^) , fi(x^) 



K " K 

J (X ) is the Jacobian evaluated at {x} , and is a parameter, 

usually 1, where \ <1 represents a Newton type of iteration and 

> 1 represents a gradient type of iteration. 

The device equations describing the behavior of a GaAlAs solar 
cell were able to fit the uniradiated data, and typical results 
are illustrated in table 6. 








Table 6. Typical results frcm device equations, 


CELL =. 35B FLU =0 

KG = 0 KP = 0 

SA = 1488265 SG = 16047.07 

D = 4.000000000E"5 DL = 0.0004 

MUP = 300 MUN = 337 3.194 

ND = 1.300000000E17 UA = 1 . OOOOOOOOOE19 

LA = 9. 000000000e“6 LG = 0.0002764919 LP = 0.00.05966502 

XP = 0.8341 Blim = 0.75 

TA = 1.000000000E~9 TG = 8 . 775881867e”10 TP = 4.594987145E 8 
DA = 0.081 DG = 87.11121221 DP = 7.747387532 
W = 1.118912755E~5 NI = 4246987.371 


SFl = 0.643F2 

= 0.4751846961 



WAVELENGTH 

MEAS RESP 

PRED RESP 

PERCENT ERROR 

0.4500 

0.1700 

0.1879 

9.5059 

0.5000 

0,5400 

0.5518 

2.1438 

0.5300 

0.5700 

0.5751 

0.8938 

0.6000 

0.5800 

0.5821 

0.3565 

0.6500 

0.5700 

0.5364 

2.8010 

0.7000 

0.5900 

0.5932 

0.5351 

0.7500 

0.6100 

0.6035 

-1.0755 

0.8000 

0.6100 

0.6188 

1.4211 

0.8500 

0.6400 

0.6400 

0.0000 

SUM OF ERRORS 

SQUARED = 107 . 

1916172 



VOPEN2 

. JSC = 14.1531103 
VOC = 0.9859651972 



i 
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When the computer was asked to minimize errors between the 

predicted and measured values for the irradiated solar cell data, 

usually one or more of the variables moved to a boundary 

limit X . or x . . , and so a separate analysis was performed 

rmax imin 

on the data by letting only one variable change at a time. 

Figures 10, 11, and 12 illustrate the typical spectral response 
curves before irradiation and after irradiating the solar cells 
with 10^® electrons/cm^ . Note that the after irradiation spectral 
response has been greatly reduced in the blue region (i.e. reduced 
more than that at the band edge) . This reduction can be partially 
explained by a reduction of the minority carrier diffusion length 
Lg in the p-GaAs region. Figures 5, 6, and 7 illustrate the 
curves which result by changing the diffusion length parameter 
and holding all the ether parameter values fixed. The effect 
that the changing of the diffusion length has on the spectral 

response at a fixed wavelength is illustrated in figure 13. 

Note that the change in the diffusion length parameter explains 

the reduced spectral response in the red region, but a change in 
at least one additional parameter will be needed to explain the 
reduced spectral response in the blue region. 

Candidates for the additional parameter change to explain 
the reduced blue response are the parameters = recombination 
velocity and = mobility. These parameter changes and their 
effects on the spectral response curve are illustrated in figures 
14, 15, and 16. Figure 16 shows the change in the spectral 
response at a fixed wavelength as the mobility changes from 3300 
cm^/volt-sec tc 5 cm^/volt-sec with the diffusion length held 
constant. This type of a reduction in the mobility is not 
realistic, and the change in the spectral response is not large 
enough to explain the overall reduced blue response. Consequently, 
the mobility parameter was ruled out. 

The recombination velocity at the p^GaAlAs , p-GaAs inter- 

face was varied, and the corresponding spectral response curves . 
at a fixed wavelength are illustrated in figures 14 and 15. 

Those curves are for various fixed values of Lg. Computed 


BEFORE IRRADIATION 





Figure 10. Spectral response curves before and after 
irradiation with 10^® electrons/cm'. . 












RECOMBINATION VELOCITY S Icm/sec) 













spectral response curves using values of S 


higher than the 


before irradiation values are also illustrated in figures 
10, 11 and 12. Note that these changes in will reduce the 

blue response and not greatly effect the spectral response in 
the red region. 

The computer analysis of these parameter changes affected 
by 1 MeV electron irradiation resulted in the paper (ref. 18) . 

A summary of the results in that paper is that the changes in 
the interface recombination velocity have been shoWn to be related 
to changes in the diffusion length produced by electron irradiation. 
These changes can be expressed 


= S 




go\ L 


(84) 


where the subscript o denotes initial values of the parameters 
before electron irradiation and S , L_ are values of the 
recoinbination velocity and diffusion length respectively after 
irradiation. 

Model for Radiation Damage 

Using the device equations developed in the previous sections, 
the computed spectral response curves can be made to agree with 
the experimental data by assuming a change in the minority carrier 
diffusion length given by (ref. 12); 


+ Kg <}. 


(85) 


where 


is the initial diffusion length (cm) , L is the final 


diffusion length (cm), 4i is the electron fluence (e/cm^) , and 
Kg is a damage coefficient. This type of change in the minority 
carrier diffusion length will account for the reduced spectral 
response in the red region of the spectrum. In addition to the 
minority carrier diffusion length change there is also a change 


. /- 


in the interface recombination velocity given by equation (84). 
This change can be related to the damage constant by using 
equation (85) 


S = S 
9 9. 




9 9r 


where is the initial value of the interface recombination 

^o 

velocity amd is the value after irradiation with a fluence 

(p. Using the above model for radiation damage together with the 
device equations yields the theoretical curves of figures 17 and 
18. These curves have the same character as the experimantal 
data obtained from irradiating GaAlAs/GaAs heteroface solar cells 
with fluences of 10^^, 10^“, 10^^ and 10^® 1 MeV electrons. 

For shallow junction cells there will also be some damage in 
the base region, and we assume 


where K is a damage constant for the base region. 

p - ^ , 

The radiation deimage model was applied to GaAlAs/GaAs solar 
cells that had been exposed to 1 MeV electron fluences of 10^^, 
lO'**, 10 and 10 electrons/cm^ . Estimates of the initial values 
of the solar cell parameters were obtained from the computer 
minimization discussed earlier. These estimates in the initial 
values were not very accurate, as can be seeri from figure 19 which 
is a numerical error analysis on the computed initial values of 
Lg and Sg. Despite the errors in the initial values an attempt 
was made to estimate a range of values for the damage coefficient 
Kg. For the cells that were tested the average damage constant 
was Kg = 7.51(10“®) with a standard deviation of ±6.93(10"®). 

This large standard deviation is due to the errors in the initial, 
fitting of the solar cell parameters. 
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UN RADIATED 


0.001 


0.0001 


Figure 




WAVELENGTH A® 

17. Theoretical spectral response curves from the 
radiation damage model and device equations . 





[mAmp/cm^ ) 



3.97 0.975 0.98 0.985, 0.99 0.995 1.0 1.05 

Vqc (volts) 

Figure 19. Error analysis on coniputed initial values of L 



The effect of the damage constant K on the spectral 

Cr 

response curves is to lower the spectral response near the band 
edge. Since we were primarily trying to model deep junction cells 
which were not the kind one would use in an actual space experiment, 
the damage coefficient Kp was kept at zero. 


Degradation of Electrical Properties 

For GaAlAs/GcLAs heteroface junction cells the quantum efficiency 
Q(X) given by equation (62) is the ratio of the photocurrent 
collected at wavelength \ to the number of photons incident 
upon the surface of the cell at wavelength \. To a first approx- 
imation we will neglect the terms occurring in 

u * ^ j ” 

equation (62). Also in equation (40) we will neglect the term 
and assiome e ' j is very small. Then the collection efficiency 
before irradiation can be approximated by: 


aLg e 

= Q^(o.,L ,x.,s ,t \ = 

^o 


^^g L 

° gp 

- L ° coshIj-L 

\ / V ^o 


(the zero subscript denoting initial value) 

The collection efficiency after irradiation with fluence is 

then 

Q*— Q*(a,L/X.,S,T) (i 

? 4> • g 1 g . .g 


where L and S after irradiation are qiven bv equations (S5) 

g g - . ' ' 

and (86). Further, since L, is related to lifetLme r by Lj 
Dr we can write 

g g 
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+ D K ({i 

g g ■ 


( 89 ) 


we make the additional assumption that the change in short-circuit 
current rates is equal to the quantum efficiency ratio before and 
after irradiation. Then for a large absorption coefficient a we 
have the approximation 



(90) 


Using equations (85) , (86) , (89) , and (90) we assign the 
nominal values K = 3.5(10""®), T * 1.22(10“®), S =10"*, 

L = 3.25(10”'*), D = 86. We can then vary the fluence d and 
% ^ 

obtain figure 20. ’ 


This figure shows that the short-circuit current ratio is 
a function of fluence and junction, depth. The curves in figure 
20 are only approximate values, but they have the same general 
shape indicated by the experimental data. 

The open-circuit voltage is related to the short-circuit 
current by relation 


V 


oc 



(91) 


where A is the perfection factor and is the dark current. 

We may write 
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is I 


V 

oc I 

‘^sc 

.^0 

■ 

+ 1 

in 

V., ~ r 

£„[. 

I 

SCO 

1 


Cn| 

L 

0 



where q = /I 

SC^ SCO 

Expanding 1 


wuefi. n - ■‘■sc'^-'-.gco’ Expanding tne numerator in equation (92) 
as a Taylor series about n = i we obtain 


I 


y = 1 +Ki(n - 1 ) + Ki(n- 1 )- + ... + K (n - i)" + ... (93) 

oco 


where Kii Kj , . . . are constants. We can also express equation 
(93) in the form 


». * + »:f^V + .. 

\ SCO/ \ SCO/ 


where a®, c\i, a;,... are consta-its. This later equation also 

suggests a quadratic term in the statistical model for V /V 

. oc oco 

as was developed earlier. 




CONCLUSIONS 

From the experimental and theoretical studies of GaAs solar 
cell behavior when exposed to 1 MeV electron f luences the following 
conclusions can be arrived at: 

1. The models for ^gc^^sco '^oc' ^oco' developed for 

predictions of degradation due to effects of 1 MeV electron f luences, 
provide excellent fits to the data and will most likely provide 
good predictions within the iexperimental region (1.5 < X . < 4.0) 
and (10*^ 2. 1 10^®). .\s coll production procedures improve the 

experimental error *,.*rll undoubtedly decrease and the derived 
prediction equations may be updated at that time, as the more 
improved data will reduce the prediction intervals and also the 
residual sum of squares of error will be reduced. 


•i ‘i « 
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2 . It must be emphasized . that the limited number of cells 

study were only test calls which were not radiation 
dened. Hence, in the use of the prediction equations to 
estunate degradation effec^s in a 

will i-fhnf. / n ertects in a space environment, the predictions 

xll tend to overestimate the actual damage. Also effects of 

nnealing are not incorporated into the prediction equations. 

tha/' ^nd experimental results indicate 

radiation hardened ceils will be these cells with small 
junction depths (in the range from 0.2 to 0.5 pm). . 

ths studies and experimental studies on damage to 

the test cells indicate that 1 „av electron fluences decrease the dif- 

isnotii Li in the n— pAAf: >-o. ^ 

ino thf. r K- ^ • region while simultaneously increas- 

ihg the recombination velocity at the GaAlAs/GaAs interface. 

sola^i t tPhta experiment to actually test GaAs 

sidered? ^ ^'"'“Pnment the following items should be con- 

(a) Goals of experiment; 

(b) Effects of shielding on electrical performance of solar 

cells in space; ^ 

r '■ 

(c) Effects of annealing on radiation damage and eventual 
electrical performance in space; 

(d) Sample size to be tested in order to achieve significant 
Statistical results; 

(e) Accurate measurement of the parameters which affect 
results (i.e., measurements of equivalent! MeV electrons, 
temperature and annealing time); and 

(f) Eff.acts Of annealing on electrical contacts and determination 

of probability of solar cell faiiivrc ano *. 

^a^lure due to external causes 

.as this probability will affect sa.mpie size on future 

. experments in space. 






FIEFERENCES 

1. Woodall, J.M.; and Hovel, H.J.: Appl. Fhys. Letters, vol. 30, 

no. 492, 1977. 

2. Draper, N.R.; and Smith, M.: Applied Regression Analysis. 

John Wiley & Sons, Inc. (N.Y.), 1966. 

3. Cochran, W.G.; and Cox, G.M. r Experimental Designs, 2nd 

Edition. John Wiley & Sons, Inc. (N.Y.), 1957. 

$ 

4. Miller, R.G. : Simultaneous Statistical Inference. McGraw- 

Hill, Inc. (N.Y.), 1966. 

5. Bancroft, T.A.: Topics in Intermediate Statistical Methods. 

The Iowa State Univ. Press (Ames , ‘ Iowa) , 1968. 

6. Netter, J.; and Wasserman, W. : Applied Linear Statistical 

Models. Irwin (Homewood, 111.), 1974. 

7. Hilberg, R. ; Teaque, M.J.; and Vette, J.I.: Comparison 

of the Trapped Electron Models AE-4- and AE-5 with AE-2 
and AE-3. NASA, Goddard Space Flight Center, NSSDC 74-13,. 

Sept. 1974. 

8. Walker, G.H.; Conway, E.J.; Annealing of GaAs Solar Cell 

Damage by Electron Irrcidiation. J. Electrochem. Soc., vol. 125, 
no. 4, Apr. 1978, pp. 676-677. 

9. Hovel, H.J.: Semiconductor and Semimetals. Vol. II: Solar 

Cells. Academic Press (N.Y.), 1975. 




■ !.g 


is ■ 

y- 

I 


■■I 












1 

'i 

r 






10. Hovel, H.J. and Woodall, J .M. : Theoretical and Experimental 

Evaluations of Ga^_j^lj^s-GaAs Solar Cells. Rec. 10th 
IEEE Phot. Spec. Conf., Palo Alto, Calif., 1973. 

11. The Solar Cell Array Design Handbook, vol. 1. Jet Prop. Lab., 

NASA, Oct. 1976. 

12. Chaffin, R.J.; Microwave Semiconductor Devices; Fundamentals 

and Radiation Effects. John Wiley & Sons (N.Y.), 1973. 

13. Hovel, H.J. and Woodall, J.M.: Optimization of Solar Cells 

for Air Mass Zero Operations and a Study of Solar Cells 
at High Temperatures. Final Report, NASA Contract NASI- 
12812, NASA-CK-145008 , 1976. 

69 . 


■i ! il 


Streetman, B.G.: Solid State Electronic Devices. Prentice 

Hall, Inc. , 1972. 

Corbett, J.W.: Electron Radiation Damage in Semiconductors 

and Metals. Academic Press (N.Y.), 1966. 

Larin, F. : Radiation Effects in Semiconductor Devices. Jo: 

Wiley & Sons, Inc. (N.Y. ), 1968. 

Cooley, W.C.; and Janda, R.J.; Handbook of Space Radiation 
Effects on Solar Cell Power Systems. NASA, SP-3003, 1963. 

Walker, G.H.; Byvik, C.E.; Conway, E.J.; Heinbockel, J.H.; 
Doviak, M.J.: Analytical and Experimental Study of 1 MeV 

Electron Irradiated GaAlAs/GaAs Heteroface Solar Cells. 
(Submitted to Electrochemical Society) 













